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Abstract — We analyze the asymptotic behavior of discrete- 
time, Markovian quantum systems with respect to a subspace 
of interest. Global asymptotic stability of subspaces is relevant 
to quantum information processing, in particular for initializing 
the system in pure states or subspace codes. We provide a linear- 
algebraic characterization of the dynamical properties leading to 
invariance and attractivity of a given quantum subspace. We 
then construct a design algorithm for discrete-time feedback 
control that allows to stabilize a target subspace, proving that 
if the control problem is feasible, then the algorithm returns an 
effective control choice. In order to prove this result, a canonical 
QR matrix decomposition is derived, and also used to establish 
the control scheme potential for the simulation of open-system 
dynamics. 

Index Terms — Quantum control, QR decomposition, invariance 
principle, quantum information. 

I. Introduction 

SINCE the pioneering intuitions of R. P. Feynmann lfT4l . 
Quantum Information (QI) has been the focus of an 
impressive research effort. Its potential has been clearly 
demonstrated, not only as a new paradigm for fundamental 
physics, but also as the key ingredient for a new generation 
of information technologies. Today the goal is to design and 
produce quantum chips, quantum memories, and quantum 
secure communication protocols fl2l . E5ll . j4|. The main 
difficulties in building effective QI processing devices are 
mainly related to scalability issues and to the disruptive action 
of the environment on the quantum correlations that embody 
the key advantage of QI. Many of these issues do not appear 
to be fundamental, and their solution is becoming mainly an 
engineering problem. 

Most of the proposed approaches to realize quantum infor- 
mation technology require the ability to perform sequences of 
a limited number of fundamental operations. Two typical key 
tasks are concerned with the preparation of states of maximal 
information ifTTl . E51 . 04ll and engineering of protected 
realization of quantum information 11351 . 11221 . Ifl8ll . IfTTl . i.e. 
the realization of information encodings that preserve the 
fragile quantum states from the action of noise. This paper 
will focus on these issues, providing a design strategy for 
engineering stable quantum subspaces. 
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In [28 1 and [ 30 1 the seminal linear-algebraic approach of 
ll22l to study noise-free subspaces has been extended to the 
general setting of noiseless subsystems |fl8l (which usually 
entails an operator-algebraic approach) and developed in two 
different directions, both concerned with the robustness of the 
encoded quantum information. The first ([28]) studies the cases 
in which the encoded information does not degenerate in the 
presence of initialization errors, the other ( [ 30 1 ) aims to ensure 
that the chosen encoding is an invariant, asymptotically stable 
set for the dynamics in presence of the noise. The latter tightly 
connects the encoding task to a set of familiar stabilization 
control problems. 

Feedback state stabilization, and in particular pure-state sta- 
bilization problems, have been tackled in the quantum domain 
under a variety of modeling and control assumptions, with 
a rapidly growing body of work dealing with the Lyapunov 
approach, see e.g. ESI, 0, ED, 02, 113, EU, ED, ED, 
[36 1, [37), and references therein. Here we embrace the 
approach of l30l , ||3T| . extending these results to Markovian 
discrete-time evolutions. 

A good review of the role of discrete-time models for 
quantum dynamics and control problems can be found in Q, 
to which we refer for a discussion of the relevant literature 
which is beyond the scope of this paper. In fact, we will 
assume from the very beginning discrete-time quantum dy- 
namics described by sequences of trace-preserving quantum 
operations in Kraus representation [19], [25 1. This assumption 
implies the Markovian character of the evolution GOl . which, 
along with a forward composition law, ensures a semigroup 
structure. We introduce the class of dynamics of interest 
and the relevant notation in Section [II] A basic analysis of 
kinematic controllability for Kraus map has been provided in 
& 

After introducing the key concepts relative to quantum 
subspaces and dynamical stability in Section III Section IV is 



devoted to the analysis of the dynamics. The results provide 
us with necessary and sufficient conditions on the dynami- 
cal model that ensure global stability of a certain quantum 
subspace. We employ a Lyapunov approach, exploiting the 
linearity of the dynamics, as well as the convex character of 
the state manifold. Lyapunov analysis of quantum discrete- 
time semigroups has been also considered, with emphasis on 
ergodicity properties, in |[6l . 

The control scheme we next consider modifies the under- 
lying dynamics of the system by indirectly measuring it, and 
applying unitary control actions, conditioned on the outcome 
of the measurement. If we average over the possible outcomes, 
we obtain a new semigroup evolution where the choice of the 
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control can be used to achieve the desired stabilization. We 
make use of the generalized measurement formalism, which 
is briefly recalled in Appendix [A] This control scheme can 
be seen as an instance of discrete-time Markovian reservoir 
engineering: the use of "noisy" dynamics to obtain a desired 
dynamical behavior has long been investigated in a variety of 
contexts, see e.g. (26), 0, OH, ll33l 



(CP) maps, that via the Kraus-Stinespring theorem |[T9l admit 
a representation of the form 



The synthesis results of Section VI include a simple charac- 
terization of the controlled dynamics that can be enacted, and 
an algorithm that builds unitary control actions stabilizing the 
desired subspace. If such controls cannot be found, it is proven 
that no choice of controls can achieve the control task for the 
same measurement. The main tools we employ come from 
the stability theory of dynamical systems, namely LaSalle's 
Invariance principle ll2D . and linear algebra, namely the QR 
matrix decomposition fl6l . We shall construct a "special form" 
of the QR decomposition: In particular, we prove that the 
upper triangular factor R can be rendered a canonical form 
with respect to the left action of the unitary matrix group. 
This result and the related discussion is presented in Section 

m 

II. Discrete-time quantum dynamical semigroups 

Let I denote the physical quantum system of interest. 
Consider the associated separable Hilbert space Hi over 
the complex field C. In what follows, we consider finite- 
dimensional quantum systems, i.e. dim(7ij) < oo. In Dirac's 
notation, vectors are represented by a ket \tp) £ Hi, and 
linear functionals by a bra, G H\ (the adjoint of Hi), 
respectively. The inner product of \ip) is then represented 
as (ip\(p). 

Let 03 (Hi) represent the set of linear bounded operators on 
Hi, $j(Hi) denoting the real subspace of hermitian operators, 
with I and O being the identity and the zero operator, 
respectively. Our (possibily uncertain) knowledge of the state 
of the quantum system is condensed in a density operator, or 
state p, with p > and Trp = 1. Density operators form a 
convex set D(Hi) C S)(Hi), with one-dimensional projectors 
corresponding to extreme points (pure states, puj = \tp)(i/j\). 
Given an X G S)(Hi), we indicate with ker(X) its ker- 
nel (0-eigenspace) and with supp(AT) := Hi ker(X) its 
range, or support. If a quantum system Q is obtained by 
composition of two subsystems Qx, Q2, the corresponding 
mathematical description is carried out in the tensor product 
space, H12 = H\®H2 1271 . observables and density operators 
being associated with Hermitian and positive-semidefinite, 
normalized operators on H\2, respectively. 

In the presence of coupling between subsystems, quantum 
measurements (see Appendix [A}, or interaction with sur- 
rounding environment, the dynamics of a quantum system 
cannot in general be described by Schrodinger's dynamics: 
The evolution is no longer unitary and reversible, and the 
formalism of open quantum systems is required |9j, 0, HI, 
ll25l . An effective tool to describe these dynamical systems, of 
fundamental interest for QI, is given by quantum operations 
||25l , |fl9l . The most general, linear and physically admissible 
evolutions which take into account interacting quantum sys- 
tems and measurements, are described by Completely Positive 



T[p]=yM kP M l 



k 



(1) 



(also known as operator-sum representation of T), where p is 
a density operator and {Mk} a family of operators such that 
the completeness relation 



E 



mIm, = 1 



(2) 



is satisfied. Under this assumption the map is then Trace- 
Preserving and Completely-Positive (TPCP), and hence maps 
density operators to density operators. We refer the reader 
to e.g. (T|, E5l . ||5l , ||9l for a detailed discussions of the 
properties of quantum operations and the physical meaning 
of the complete-positivity property. 

One can then consider the discrete-time dynamical semi- 
group, acting on D(Hi), induced by iteration of a given TPCP 
map. The resulting discrete-time quantum system is described 
by 

p(t + 1) = T[p(t)} = £ M k p(t)Ml (3) 

k 

Given the initial conditions p(0) for the system, we can then 
write 

p(t)=T t [p(Q)} t = l,2,... 

where T'[-] indicates t applications of the TPCP map T[-j. 
Hence, the evolution obeys a forward composition law and, in 
the spirit of HI, is called a Discrete-time Quantum Dynamical 
Semigroup (DQDS). Notice that while the dynamic map is 
linear, the "state space" D(Hi) is a convex, compact subset 
of the cone of the positive elements in S)(Hi). 

While a TPCP maps can indeed represent general dynamics, 
assuming dynamics of the form Q, with M^'s that do not 
depend on the past states, is equivalent to assume Markovian 
dynamics (see [20 1 for a discussion of Markovian properties 
for quantum evolutions). From a probabilistic viewpoint, if 
density operators play the role of probability distributions, 
TPCP maps are the analogue of transition operators for clas- 
sical Markov chains. 

III. Quantum subspaces, invariance and 

ATTRACTIVITY 

In this section we recall some definitions of quantum 
subspaces invariance and attractivity. We follow the subsystem 
approach of [30], [31], focusing on the case of subspaces. This 
is motivated by the fact that the general subsystem case is 
derived in the continuous-time case as a specialization with 
some additional constraints, and that for many applications of 
interest for the present work, namely pure-state preparation 
and engineering of protected quantum information, the sub- 
space case is enough, as it is suggested by the results in [ 30 1 . 

Definition 1 (Quantum subspace): A quantum subspace S 
of a system I with associated Hilbert space Hi is a quantum 
system whose Hilbert space is a subspace Hs of Hi, 



Hi = Hs © Hr, 



(4) 
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for some remainder space Hr. The set of linear operators on 
S, B(Hs), is isomorphic to the algebra on Hi with elements 
of the form Xi = X s © O r . 

Let n = dim(7i/), m = dim(Hs), an d r — dim(7i/?), 
and let {\4>)j}^i, {\4>)k'}k=i denote orthonormal bases for 
Hs an d Hr, respectively. Decomposition Q is then naturally 
associated with the following basis for "Hf. 

{lw)} = {|0)f}J=iU{|0)I?}U- 

This basis induces a block structure for matrices representing 
operators acting on Hr. 



X = 





X P ' 


. X Q 


Xr 



In the rest of the paper the subscripts S, P, Q and R will 
follow this convention. Let moreover II5 and Hr be the pro- 
jection operators over the subspaces Hs and Hr, respectively. 
The following definitions are independent of the choices of 

Mf}r =1 . {\muv 

Definition 2 (State initialization): The system X with state 
p 6 T>(Hi) is initialized in S with state ps € T>(Hs) if p is 
of the form 

r II 1 

(5) 



Ps 


" 









p 

We will denote with Js(Hi) the set of states of the form |5]) 
for some ps £ V(Hs)- 

Definition 3 (Invariance): Let X evolve under iterations of 
a TPCP map. The subsystem S supported on the subspace 
Hs of Hi is invariant if the evolution of any initialized p € 
Js(Hj) obeys 





Pit) 











G Js(Hi) 



Vi > 0, and with T s being a TPCP map on H S - 

Definition 4 (Attractivity): Let X evolve under iterations of 
a TPCP map T. The subsystem S supported on the subspace 
Hs of Hi is attractive if Vp £ T>(Hi) we have: 

Km ||T*(p)-nsT*[p]n s || =0. 

t — >oo 

Definition 5 (Global asymptotic stability): Let X evolve 
under iterations of a TPCP map T . The subsystem S supported 
on the subspace Hs of Hi is Globally Asymptotically Stable 
(GAS) if it is invariant and attractive. 

IV. Analysis Results 

This section is devoted to the derivation of necessary and 
sufficient conditions on the form of the TPCP map T for a 
given quantum subspace S to be GAS. We start by focusing 
on the invariance property. 

A. Invariance of 3 s^Hi) 

The following proposition gives a sufficient and necessary 
condition on T such that JsiHi) is invariant. 

Proposition 1: Let the TPCP transformation T be described 
by the Kraus map ([TJ. Let the matrices M k be expressed in 
their block form 

M \M kiS M kt p 
k [M ktQ M ktR 



according to the state space decomposition Q. Then the set 
Js{Hi) is invariant if and only if 



Vfc. 



(6) 



Proof: Verifying Definition [3] is equivalent to verifying 
that there exists a TPCP map Ts such that 



T 



Ps 




Ts{ps) 




(7) 



for all ps in D(Hs)- By exploiting the block form of the M k 
matrices in ([TJ given by the decomposition we have 



T 



Ps 




k 

= E 



Ps 


k,S 



k L 

E 

k 






M , 

M, 



Ml 



Mk, 

M k,sPsM. 



k,R 
t 





Ps 


0" 















Ks 



rt 



k,Q 



k,R 



k.S 



M. 



k,QPsM kS 



M k.sPsM KQ 
M k, Q PsMl Q 



(8) 

Sufficiency of |6]) to have invariance of JsiHi) is trivial. 
Necessity is given by the fact that the lower right blocks 
M k QpsMl q are positive semi-definite for all fc's, and there- 
fore, for to hold, it has to be M k Q psM\ q — Vfc. For 
M k QpsM k q to be zero for any state ps € V{Hs), it has then 
to be M k q = 0. Equation Q then implies that the completely- 
positive transformation 

T s[ps] 



Y, M k,sPsMls 



is also trace preserving. 



B. Global asymptotic stability of \Js(Hi) 

The main tool we are going to use in deriving a characteriza- 
tion of TPCP maps that render a certain Hs GAS, is LaSalle's 
invariance principle, which we recall here in its discrete time 
form EJJ. 

Theorem 1 (La Salle's theorem for discrete-time systems): 
Consider a discrete-time system 

x{t+l)=T[x(t)] 

Suppose V is a C 1 function of x € W L , bounded below and 
satisfying 



AV(x) = V(T[x]) - V(x) < 0, Vx 



(9) 



i.e. V(x) is non-increasing along forward trajectories of the 
plant dynamics. Then any bounded trajectory converges to 
the largest invariant subset W contained in the locus E = 
{x\AV(x) = 0}. 

Being any TPCP map a map from the compact set of 
density operators to itself, any trajectory is bounded. Let's 
then consider the function 



V{p) = Tr(n fl p) > 0. 



(10) 



The function V(p) is C 1 and bounded from below, and it is 
a natural candidate for a Lyapunov function for the system. 
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In fact, it represents the probability of the event IT^ (see 
Appendix |Aj, that is, the probability that the system is found 
in Hr after the measurement. 

Lemma 1: Let T be the generator of a DQDS, and as- 
sume a given quantum subsystems S to be invariant. Then 
V(p) = Tr (Ujip) (Hr being the remainder space) satisfies 
the hypothesis |9]l of Theorem [T] 

Proof: The variation of V(p) along forward trajectories 
of the system ([3) is 



AV(p) = Tr(U R T[p])-Tr(U RP ) 



Tr 



P 



(11) 



Notice that Tr(J2 k M kP M l ~ p) = °> and that V (p) = for 
all p's that have support in Hs- Let us express M.pMt —p 
in its block form, using the fact that Mq — by assuming 
invariance of Js(Hi). We get 

MkpMl -p = 



E 



M, 



k,S 





M 



fe.p 



k,R 



PS 
P P 



Pp 
PR. 



M 



k,R 



E 



A h,sPsMl s +M k pp^Mt s +M k S p P Ml p +M k p p R Ml r 
M k,RPP M ls+ M k,RPR M lp 



,, S PP M t,R + M k,pPR M l,R 
M k,R.PR M l,R 



ft 





" PS 


pp ' 










PR 



(12) 



Therefore 



AV(p) = Tr 



Tr 



n fl {TMkpMl-pj 

M k,RPR M t, R ~ PR 



(13) 



so that in order to get AV < the map Tr\pr\ := 
2~2k Mk,RPnMl R has to be trace non-increasing. 

Note that this condition is automatically verified, once T is 
a TPCP map. Indeed, consider the application of T on a state 
p which has support on Hr. According to the block from in 
(fT2|i we have that the total trace of T[p] is 



Tr (T[p]) = Tr ^ M^rM^ +Tr ^ M^rM^ 

Therefore, as both the terms are positive, being pn > 0, and 
as T is TP, we have for any pR 6 V(rin) 



Tr 



M k,RpR M lR 



<Tr(T[p]) = Tr(p R ) 



and thus Tr is trace non-increasing. ■ 
This leaves us with determining when Js(H.i) contains 
the largest invariant set in E. We shall derive conditions 
that ensure that no other invariant set W exists in E = 
{p\AV(p) = 0} such that J s (Hi) C W. We start by giving 
some preliminary results. 



Lemma 2: Let T be a TPCP transformation described by 
the Kraus map ([TJ. Consider an orthogonal subspace decom- 
position Hs^BH-r- Then the set jR^Jij) contains an invariant 
subset if and only if it contains an invariant state. 

Proof: The "if" part is trivial. On the other hand, Jr^Hj) 
is convex and compact, hence if it contains an invariant subset 
W it also contains the closure of its convex hull, call it W. 
The map T is linear and continuous, so the convex hull of 
an invariant subset is invariant, and so is its closure. Hence, 
by Brouwer's fixed point theorem l40l it admits a fixed point 
peWCJaCHj). ' m 

Lemma 3: Let W be an invariant subset of T>(Hi) for the 
TPCP transformation T, and let 



H 



= supp(T40 = (J supp(p). 

pew 



Then Jw("Hi) is invariant. 

Proof: The proof follows the one for the continuous-time 
case in BTl Lemma 8]. Let W be the convex hull of W. By 
linearity of dynamics it is easy to show that W is invariant 
too. Furthermore, from the definition of W, there exists a 
p G W such that supp(/5) = supp(V^) = Hw- Consider 
the decomposition Hi = Hw ® H w , and the corresponding 
matrix partitioning 



X 



Xw 
Xm 



With respect to this partition, pw is full rank while pl,m.n 
are zero blocks. The state p is then mapped by T according to 
(|8]i and therefore, as pw is full rank, it has to be Mk,Q — for 
all k's. Comparing it with the conditions given in proposition 
[T] we then infer that JwO~Li) is invariant. ■ 
Proposition 2: Consider p £ Jr(Hi) and evolving under 
the TPCP transformation T described by the Kraus map Q. 
Let the matrices M. be expressed in the block form 



M, 



Mk.P 



M k,S 
. M k,R_ 



according to the state space decomposition Hs © Hr, with 
Js(Hi) invariant. Then p e E = {p e V{Hi)\AV{p) = 0}, 
where V(p) is defined by ( fT0] >, if and only if its pr block 
satisfies 



supp(pfl) C Qker ( 



M, 



k 



Proof: By direct computation, see ( [T2| i, we have 

^ M oM ] = V M k,PP RM tp M k,pPR M t,R 

k V m krPr m Ip m kr.prM{r 



(14) 



as p has support on Hr alone. Note that as V(p) = 1, 
AV(p) = is equivalent to T[p] having support on Hr. 
Given the form of the upper-left block of ( fl4| ), this is true if 
and only if supp(pR) C f] k ker (^M k . ■ 
Proposition [2] allows then to state the following key char- 
acterization of global, asymptotical stability of Js(Hi). 
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Theorem 2: Let the TPCP transformation T be described 
by the Rraus map ([TJ. Consider an orthogonal subset decom- 
position Hs © Hr, with Jsijii) invariant. Let the matrices 
M k be expressed in their block form 



AI, 



k,S 




M l 

M, 



k,P 



k,R 



according to the same state space decomposition. Then the set 
Jsi'Hi) is GAS if and only if there are no invariant states 
with support on 

f|ker (M KP ). 

k 

Proof: Necessity is immediate: if there was an invariant 
state with support on f] k ker I M k p J , It would have non triv- 
ial support on Hr, and therefore Jig could not be attractive. 
In order to prove the other implication, consider LaSalle's 
theorem. By hypotesis, JsiHi) is invariant and is contained 
in E (see Proposition [2|, therefore it is contained in the largest 
invariant set W in the zero-difference locus E. 

Let us suppose that J s {Hi) C W, but J s {Hi) j= W. That 
is, there exists a set W C E which is invariant and strictly 
contains Js^Hi)- Therefore its support has to be 

Ti-w — Hs © Ti.Ri 

with Hr> subspace of Hr, and by Lemma [3] Jwi'rii) must 
be invariant too. Consider then a state p which belongs to 
Jw(T~Li), with non trivial support on Hr', and define 





PR 



r Tr{U R ,p) 

which has support on Hr> only. By construction, p is in 
Jw(7~Li)> an d therefore its trajectory is contained in Jw(7~Li)- 
It is also in E, that is AV(p) — 0. As we have V(p) = 1, 
then its evolution must be also remain in TCr' C Hr at any 
time. Therefore an invariant set with support on Hr, exists. By 
reversing the implication, this means that if does not exist an 
invariant set with support on Hr, then Js{Hi) is the largest 
invariant set in E. Furthermore, Proposition [2] indicates that if 
there is an invariant set in E with support on Hr, its support 
must actually be contained in 



n 



ker M 



L k,P 



Therefore, if no such subset exists, we have attractivity of 
■Js{T~Li) by LaSalle's theorem. By Lemma|2] the existence of 
an invariant set is equivalent to the existence of an invariant 
state with support on f] k ker (^M k p^j . ■ 
Given the usual decomposition Hi = Hs © Hr, let us 
further decompose Hr, in 

H RI = H R p| ker fad h p ) and Hr» = f] ker (M fcjP ) 

k k 

and consider the operation elements M k in a basis induced by 
the decomposition Hi = Hs © Hr> Hr»,: 





' M k,S 














M k , R2 







M k,R3 





Density operators p which have support on the bottom right 
block clearly belong to Jr"(Hi). Sufficient, although not 
necessary, condition to be sure that no invariant sets have 
support on that subspace is that 



f>er( 



M, 



k.R2 



= {0} and 



This way, the states that have support on P| fc ker^Af fep ^ 
will be mapped into states which has non-trivial support on 



Dfe ker ( M, 



k,P 



and therefore no invariant set will exist in 



f] k ker ( M k p ) . This intuition will be further developed in the 



Section VI where a control design tool capable of achieving 



attractivity of a given subspace is obtained. 

V. A CANONICAL MATRIX FORM BASED ON THE QR 
DECOMPOSITION 

A. On the uniqueness of QR decomposition 

Definition 6 (QR decomposition j|76l/): A QR decomposi- 
tion of a complex-valued square matrix A is a decomposition 
of A as 

A = QR, 

where Q is an orthogonal matrix (meaning that Q^Q = I ) 
and R is an upper triangular matrix. 

We investigate here the uniqueness of the QR decomposition 
of a complex-valued matrix A, both in the case in which A 
is non-singular and in the case in which it is singular. While 
the real-matrix case is well known (see e.g. Ifl6l ). a little extra 
care is needed in the complex case. Let's first analyze the case 
of A non-singular. We have the following result. 

Lemma 4: Let 



A = Q 1 R 1 



Q2R2 



be two QR decompositions of the same non-singular square 
matrix A. Then R 1 and R 2 only differ for the phase of their 
rows, that is 



R l = $i?2 . 



j'i'i 



Qi = Q2®~ 



c,3<kn 



where $ = diag (e 

Proof: Using the fact that A and its factors are non 
singular, we have 

q\q 1 = i? 2 i?r l = $ 

where $ is upper triangular because it is the product of 
two upper triangular matrices. The matrix $ must also be 
orthonormal, because it is the product of two orthonormal 
matrices. Therefore, starting from the first column of we 
have that |3>u| = 1. where <i> i:) is the element of $ in position 

We proceed by induction on the column index j. Assume 
that all the columns Q k with k < j satisfy — for any 
I ^ k. In order for 5>j to satisfy & k <&j — Vfc < j, it must 
be 



1.7 



0. 
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Moreover, as $ is upper triangular, we must also have 
$j+i j = ■ ■ ■ = j = 0. Therefore, by othonormality of 



$, it has to be <& 



33 I 



1. 



In the case in which A is singular, on the other hand, the 
QR decomposition is not just unique up to a phase of the rows 
of R. 

Example 1. Consider the following matrix: 



M = 



Since it is already upper triangular, a valid QR decomposition 
is given by Q = I, R = M. On the other hand, we can 
consider 



Q 



l 

7^ 



i 

-i 



R 



V2 




which clearly also give QR — M, with Q'Q = I. 

However, introducing some conditions on the QR decom- 
position, it is possible to obtain a canonical form for the QR 
decomposition in a sense that will be explained later in this 
section. A useful lemma in this sense is the following. 

Lemma 5: Consider a QR decomposition of a square matrix 
A of dimension n, and an index j in [l,n], such that 



nj = Vj < j,Vi > pj 



(15) 



where pj is the rank of the first j columns of A. Let and 
<7i, be the i-th column of A and Q respectively. Then 

<ai,... ) a j > = <qi > ...,q Pj > Vj = 1, . . . , j. 

Proof: Consider the expression for the j-th column of A 

a 3 = Qrj. 

Pj elements of Vj are zeros, 



By the hypothesis, the last n 
hence it results 



dj e < qi, 
and therefore 



,q P] > Vj = i,...,j 



< a-i 



,a 3 >C< q x ,...,q Pj 



> 



Vj = 1, 



As the rank of the first j columns is pj, which is also the 
dimension of < q%, . . . , q p . >, equality of the two subspaces 
holds. ■ 
We next show as a particular choice of the QR decomposi- 
tion, suggested by lemma [5] gives a canonical form on C nx ™ 
with respect to left-multiplication for elements of the unitary 
matrix group U(n). We construct the QR decomposition 
through the Gram-Schmidt orthonormalization process, fixing 
the degrees of freedom of the upper-triangular factor R and 
verifying that the resulting decomposition satifies the hypotesis 
of lemma [5] for j = n. 

B. Construction of the QR decomposition by orthonormaliza- 
tion 

Theorem 3: Given any (complex) square matrix A of di- 
mension n, it is possible to derive a QR decomposition 
A = QR such that hypoteses of Lemma [5] are satisfied, and 



such that the first nonzero element of each row of R is real 
and positive. 

Proof: We explicitly construct the QR decomposition of 
A column by column. We denote by A, Q, R the matrices, 
with a,i,qi,ri their i-th columns and with aij,qij,ri,j their 
elements, respectively. Let us start from the first non zero 



'7i 



n. 



»'2.. 



(16) 

Also fix rj = for all j < iq. 

The next columns of Q,R are constructed by an iterative 
procedure. Define p^\ as the rank of the first i — 1 columns 
of A. We can assume (by induction) to have the first pi-\ 
columns of Q and the first i — 1 columns of R constructed in 
such a way that rkj = for k > pj and j < i — 1. 

Consider the next column of A, cii. Assume that is 
linearly dependent with the previous columns of A, that is 
pi = Pi-i. Since lemma [5] applies, can be written as 



i-l 



-1 



Pi 



^ a 3 a -j = 2^ a i 2^ T ^ qi 

j = l 1 = 1 1=1 



and therefore, being a linear combination of the columns 
{qi, . . . ,q Pi l }, the elements of Ti are defined as 



for t= 1,. 



,Pi 



On the other hand, if the column ai is linearly independent 
from the previous columns of A, then the rank pi = pi-i + 1. 
As before, the first p^i coefficients of r t must be defined as 



/ , , = qjai , for £ = 1, . . . ,pi - 1, 
and let set r^i = for I = pi+1, . ■ ■ , n. Let us also introduce 



Pi 

E 

1=1 



rtjqt ^ 



and define 



Tpi,i — || 



(17) 



(18) 



It is immediate to verify that the obtained q Pi is orthonormal 
to the columns gi, . . . , q Pi -i, and that a; = Qr Pi . 

After iterating until the last column of R is defined, we 
are left to choose the remaining columns of Q so that 
the set {qi, ■ ■ ■ ,q n } is an orthonormal basis for C" xn . By 
construction, A = QR. ■ 



C. R is a canonical form 

Let Q be a group acting on C nx ". Let A,B e C nxn . If 
there exists a g S Q such that g(A) = B, we say that A and 
B are -equivalent, and we write A^g B. 

Definition 7: A canonical form with respect to Q is a func- 
tion T : C nxn -> C" x ™ such that for every A,B e C nxn : 

i. T{A) ~ Q A; 

ii. T{A) = T(B) if and only if A ~g B. 
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Let us consider the unitary matrix group U(n) C C nxra and 
consider its action on C nxn through left-multiplication, that 

is, for any U G U{n), M € C nx ": 

U(M) = UM. 

We are now ready to prove the following. 

Theorem 4: Define J 7 (A) = R, with R the upper-triangular 
matrices obtained by the procedure described in the proof of 
theorem [3] Then J 7 is a canonical form with respect to U(n) 
(and its action on C nxn by left multiplication). 

Proof: By construction A = QR, with unitary Q, so 

H A ) ~W(n) A If A , B G C " X " ^ such that -T 7 ^) = 
^"(JS) = i?, thus A = QR and _B = Vi? for some 

Q, V e U{n), and hence A = QV^B. 

On the other hand, if A — UB, U € U{n), we have to 

prove that the upper-triangular matrix in the canonical QR 

decompositions A = QR^ and B = VR^ is the same. If 

the first non zero column of B is 6j , then the first column 

nonzero column of A is, being U unitary, a, = Ubi . One 

then finds from ( fTo*) 

Hence the first irj columns of and R^ are identical. We 



9i 



(A) 



then proceed by induction. Assume that ri- 
f/wj for j = 1, i — 1. If the column a, is linearly dependent 
from the previous i — 1 so it must be 6, . The elements of r 
are defined as 

r i A } = QkOi = qlurfai = vhi = r [ k B } , for k = 1, 



On the other hand, if the column a% is linearly independent 
from the previous columns of A, then the rank pi — + 1. 
As before, the first Pi — 1 coefficients of are defined as 



k.i 



q\a.i 



and = r K 

as before 



= for k 



r k,i 1 



for k = 1, 



■ A' 



, n. Let us consider 



fc=i 



By using the equivalent definition and the inductive hypothesis 
it follows that b L = Wdi and 



\U*Oi\ 



Hence r\ A ^ — r^ B \ and by induction R^ = R( B \ 



in an open system, discrete-time dynamics described by a 
given Kraus map, with associated Kraus operators {Mk}- This 
can be realized, for example, when the system is coupled 
to an auxiliary measurement apparatus, it is manipulated 
coherently, and then a projective measurement is performed 
on the auxiliary system (see appendix [Aj. Suppose moreover 
that we are allowed to unitarily control the state of the system, 
i.e. Pcontroiied = U pW, U € U{TLi). We shall assume that 
the control is fast with respect to the measurement time scale, 
or the measurement and the control acts in distinct time slots. 

We can then use the generalized measurement outcome k 
to condition the control choice, that is, a certain coherent 
transformation Uk is applied after the fc-th output is recorded. 
The measurement-control loop is then iterated: If we average 
over the measurement results at each step, this yields a 
different TPCP map, which dependens on the design of the 
set of unitary controls {Uk} and describes the evolution of 
the state immediately after each application of the controls: 



p(t + l) = J2u k M k p{t)Mlul 



Figure [T] depicts the feedback control loop (before the averag- 
ing). 



Gen. Measurement 



Unitary Control 



Pi 



{M k } 



Pt\k 



M kPt Ml 



{Ut} 



trace (M kPt M}\ 
k 



Pt+i = U k pt\ k Ul 



-1 



Classical Information 
Fig. 1. A measurement-dependent unitary control scheme. 

In this section we will tackle the problem of characterizing 
the set of open loop dynamics that can be engineered through 
this feedback setup by designing the set {Uk} with fixed 
measurement operator {Mk}. Moreover, on the basis of the 
analysis results of section [IV] we will derive an algorithm that 
allows to design the set of unitary controls such that global 
asymptotical stability of a given subspace is achieved. Notice 
that the desired result is going to be achieved for the averaged 
time-evolution that describes the system immediately after the 
control step. It is easy to show that if a certain subspace is 
GAS for the averaged dynamics, it must be so also for the 
conditional ones. 



VI. Engineering attractive subspaces via 

CLOSED-LOOP CONTROL 
A. The controlled dynamics 

In this section we deal with the problem of stabilization of 
a given quantum subspace by discrete-time measurements and 
unitary control. The control scheme we employ follows the 
ideas of ED, l29l . 

Suppose that a generalized measurement operation can be 
performed on the system at times t = 1,2,..., resulting 



B. Simulating generalized measurements 

A first straightforward application of the canonical form we 
derived in the previous section is the following. Assume we are 
able to perform a generalized measurement, with associated 
operators {Mk}u n —i, and we would like to actually implement 
a different measurement with associated operators {Nk}™^, 
by using the unitary control loop as above. Notice that the 
control scheme we considered allows to modify only the 
conditioned states, not the probability of the outcomes, since 
trace(M J jM fey o) = trace [Mlu\U k M k p). The following holds: 
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Proposition 3: A measurement with associated operators 
{N k }™ =1 can be simulated by a certain choice of unitary 
controls from a measurement {M k }™ =1 , if and only if there 
exist a reordering j(k) of the first m integers such that: 

where T returns the canonical R factor of the argument, as 
described in the section [V] 

Proof: Let us first assume that for a given reordering j(k) 
it holds !F(N k ) = ^-"(MjVfc-j) = R k . Therefore the canonical 
QR decomposition of N k and Mjr k \ gives 



N k — U k R k 



Mj(k) — Vj(k)Rk- 



Let then U k Vj^ be the unitary control associated with the 
measurement outcome k. We have 



%losed loop [p] = J2 UkV m M 3(k)P M j( k ) V 3 



m u k 



k k 

and therefore it simulates the measurement associated with the 
operators {iV fc }™ =1 . 

On the other hand, suppose that there exists a set of unitary 
controls {Qk}™ =1 and there is a reordering j(k) of the first 
m integers such that 

Qj(k)Mj(k) = N k . 

According to Theorem |4j T is a canonical form with respect 
to U(n) and its action on C nxn by left multiplication, and 
therefore if Q m M m = N k then T{M ][k) ) = T{N k ). ■ 

C. Global asymptotic stabilization of a quantum subspace 

Suppose that the operators {M k } are given, corresponding 
to a measurement that is performed on the quantum system, 
with corresponding outcomes {k}. We are then looking for a 
set of unitary transformations {U k } such that, once they are 
applied to the system, the resulting transformation 



T[p] = Y J U k M kP Mlul 



R 



makes a given subsystem S GAS. Let us introduce a prelim- 
inary result. 

Lemma 6: Let R be the upper triangular factor of a canon- 
ical QR decomposition in the form 

Rs Rp 
R R 

(according to the block structure induced by Q) and suppose 
Rp = 0. Consider the matrix N obtained my left multiplying 
R by a unitary matrix V: 



N = VR = 



Then N Q = implies N P = 0. 

Proof: Let first consider the case in which Rs has full 
rank. Let r x m be the dimension of Vq, and m x m be the 
dimension of Rs- Since it must be Nq — VqRs — and R$ 
is full rank, we have Vq = 0. 



'Vs Vp 




'Rs " 




~N S N P 


Vq Vr 




R R 




Nq Nr_ 



As V is unitary, its column must be orthonormal. Being 
Vq — 0, Vs must be itself an orthonormal block in order 
to have orthonormality of the first m columns of V. It then 
follows that Vp = 0, because any j-th column, j > m, must 
be orthonormal to all the first m columns. It then follows that 



Nr 



V p Rr 



0. 



Let us now consider the other case, in which Rs is singular. 
This implies that p m < m (p m being the rank of the first m 
columns of R): Therefore, as R is a triangular factor of a 
canonical QR decomposition, the element R Pm +i. m = 0. 

Now, by construction of the canonical QR decomposition, 
if there were non-zero columns of index j > m, one of them 
would have a non-zero element on the row of index p m + 1. 
By recalling that Rp = 0, we have that R Pm +ij = 0, Vj € 
[m + l,7n + r]. Therefore, all the last r columns are zero- 
vectors, and in particular Rr = 0. It then follows that 

N P = VrRr = 0. 



This result will be instrumental in proving the main theorem 
of the section, which provides us with an iterative control 
design procedure that renders the desired subspace asymptot- 
ically stable whenever it is possible. 

Theorem 5: Consider a subspace orthogonal decomposition 
Hi = TLs © T^-R and a given generalized measurement 
associated to Kraus operators {M k }. If asymptotic stability of 
a subspace S can be achieved by any measurement-dependent 
unitary control {U k }, it can be achieved by building U k using 
the iterative algorithm below. 



Control design algorithm 



Let {\<f>)f}f =1 , {\4>)k}l=i denote orthonor- 
mal bases for Hs and Hr, and represent each 
M k as a matrix with respect to the basis 
{\<P)j IHMfe } r k=v Compute a QR decomposi- 
tion M k = Q k R k with canonical R k for each k. Call 
?4 0) = Hr, = U s , U [ k ] = Q{ and rename 



the matrix blocks Rc \ = Rs 



Rp\ = Rp, k and 



K R,k - K R,k- 

If Rp. k = Vfc, then the problem is not feasible 
and a unitary control law cannot be found. Otherwise 
define V (0) = /, Z (0) = I, and consider the 
following iterative procedure, starting from i = 0: 



1) Define H 



= n 



if n 



(i+l) 



k ^Rp, k 



{0} then the iteration is success- 



fully completed. Go to step 8). 
If H^ +1} C ftW, define 

n% +f) and = /. 



e 



(t+1) _ 

R 



H% (i.e. R { p ] k = Vfc) then, if 



C jS of 



if 

dim(W^) > dim{Hf): 

a) Choose a subspace Hs +1 ^ 

(») 



the same dimension of TL\ . (Re)-define 



n 



(<+l) 



H 



(i) 



G^ +1) . 



9 



b) Let T-6j) = 0*~q Tig\ Consti-uct a unitary 
matrix Y with the following block form, 
according to a Hilbert space decomposition 



<1V 

ill — rLrp 



H 



H 



(i+l) 



Y 



7 








' 





1/^27 










1/V2I 


-1/V27 














7 



(20) 



If instead dim(H^) < dim(W^) 



a) Choose a subspace Tig +1 



same dimension of 7Y 



b) Let 

( 



n 



(?) 



C W^ 3 of the 



Construct a unitary 
matrix Y with the following block form, 
according to a Hilbert space decomposition 



Til = Tin 



(i+l) 



7^ 





" I 











1/V2I 


1/72/ 




_ 


1/V2I 


-1/V2/ 



(21) 

c) Define Z< 4+1 ) = zMy('+!) and go to step 
8). 

2) Define = Z { ^Y {i+1 \ 



3) Rewrite fl^ = W {i+ ^ R ( ^ k W {l+1 ^ in a 



basis according to the Tip^ = Ti^ +1 ^ i 



(i+l) 



decomposition. 
4) Compute the canonical QR decomposition of 



R 



oi: +1) ^ +1) . 



Compute the matrix 



blocks Rp^KR^ of 77^ +1) , again accord- 



ing to the decomposition Ttjj = Ti 



/(i+l) 



5) Define 

7 ' 



[/('+!) 





w (i+1)t (<9i i+1) ) f 



6) Define 



7 











7) Increment the counter and go back to step 1). 

8) Return the unitary controls Uk = 



Proof: Let us first consider the case in which the algo- 
rithm stops before the iterations. This happens if for every k 
we have Rp : k — 0. Remember that each 77fc has been put 
in canonical form, so Lemma [6] applies: This means that any 
control choice that ensures invariance of the desired subspace, 
that is Nk = UkRk with Nq^ — 0, makes also Jp(Ti[) 
invariant, since Npk — 0. Hence an invariant state with 
support on Tip, always exists. This, via Theorem [2] precludes 
the existence of a control choice that renders Js{Tii) GAS. 

If the algorithm does not stop, then at each step of the 

/(0 



hence the algorithm is completed in at most n steps. If the 
algorithm is successfully completed at a certain iteration j, we 
have built unitary controls {U^ } and a unitary such that 
the controlled quantum operation element, under the change 
of basis V^', is of the form: 



V^UkMkV^ 



Z U) 



n S,k 


n P,k 














n S,k 


















n P,k 














M 

H S,k 


5(3) 

H p,k 














K R.k 



(22) 



where the block structure is consistent with the decompo- 
sition ©1=0 (where to simplify the notation we set 



n 



O'+i) _ 



= Ti'ii')- Let 77fc be the block matrix above and 



consider its upper-triangular part. The rows have the form 



I?' 



(0 
P,k 











because at each step of the iteration 



we choose a basis according to the decomposition 



n 



(i+l) 



Ti 



(i+l) 



where Tip +1 ^ C f], ker Rp' k , hence ob 



taming R p ' k 



R 



(i) 
P.k 











It is easy to verify 



that the subsequent unitary transformations have no effects on 
the blocks Rp^ k . 

The upper-triangular form of each Rk and the form of Z^>> 
and both block-diagonal with respect to the orthogonal 
decomposition Tis © Tip, ensure invariance of Tis- 

By construction, for all i = 0, . . . , j, either f] k ker R^p k = 
{0} and Y® = 7, or R$ k = for all k and differs from 
the identity matrix and has the form ( |20] > or ( pi) . 

Let us prove that no invariant state can have support on 

7-4-1 ( -pi 

©i=i by induction. First consider a state with support 



on Ti 



O'+i) 



Ti^ alone: 



" 








PR _ 



If akeri?^ 



= {0}, then p is mapped by Y^k R k *<- k 
into a state p' with non-trivial support on Tig . Being in 
this case = 7, is block-diagonal with respect 

to the considered decomposition and we cannot thus get 

ZWpfzVH = p, for any p in D(7Y^). 



On the other hand, if 77 



(j) 



Vfc, then yCi) contains 



off-diagonal full-rank blocks and maps the state 





P 



l^k K R,kPR U R,k 



into a state with non-trivial support on Tig . The subsequent 
application of Z^^ 1 ^ will then map the state into a state with 
nontrivial support on Q^j^s', and therefore p cannot be 
invariant. 

Let us now proceed with the inductive step, with m as 
the induction index. Assume that no invariant state can have 



iteration the dimension of Tip' is reduced by at at least 1, support on @i=j+i— m Tig 1 alone (induction hypothesis), and 



to 



1 + 1 

consider the subspace (£)^ =7 _ 



Ji g ■ By the induction hy- 



pothesis if there were an invariant state with support on this 
subspace, it would be in the form 





' 








p = 





PS 


pp 




_ 


A 


PR 



with ps ^ having support on Jig . 



Let us rewrite 



g(J) — ^(3-m-l)y(j-ra) 



trivial support on Ji^j r " 1 



where Z (j _ m+1) = yO'-m+i) ..... yCO. 

Again, we have two cases. If f] k ker Rpi™ ^ = {0}, then 
yu'-ro.) = / and p is mapped by 77fc into a state with non 

The subsequent application of 
Z{j-m+i) ar >d of y( J ~ m ) does not affect this, and because of 
ZU - " 1-1 ), the first complete iteration will map p into a state 
with non trivial support on 1 Jig . Therefore p cannot 

be invariant. 

On the other hand, if R ( ^ k ' n ~ 1] = Vfc, then F^'~ TO ^ has 
the form (|20|> and the closed loop evolution of p is 



P = 



E 

k 



Z (3~ 



t)y(i-ro) ^ 



(j-m+l)RkpR k Z l 



G-m+1) 



y(i-"i)t^(i-m-i)t 



If /5fc has support on 0^=j + i_ OT W<j for all k, then p' will 
have the same support, and therefore p is not invariant. If 
instead pk has non trivial support on Jig for some k, then 
because of the subsequent application of ^U _m_1 )y(J _m ), p 
will have non trivial support on ©^ZT 1-1 He , and again p is 
not invariant. 

When the induction process reaches m = j — 1, then it states 



that no invariant states are supported on 77. 



(1) 



0) 



77^-', and therefore according to Theorem|2|global asymptotic 
stability of the subspace S is achieved. 

The algorithm is clearly constructive. We then get the 
following: 

Corollary 1: A certain subspace 77s can be made GAS 
if and only if the Rpk blocks of the canonical R-factors, 
computed with respect to the decomposition 77/ = 77s ®JiR, 
are not all zero. 



VII. A Toy Problem 

We consider in this example a two-qubit system, defined on 
a Hilbert space 77/ ~ C 2 <g)(C 2 . Consider the task of stabilizing 
the maximally entangled state 



Pd = \ (100) + |11)) «00| + (11|) 



(23) 



which has the following representation in the computational 
basis C = {\ab) = |o) ® |6)|a,6 = 0, 1}: 



Pd = 



10 1 





10 1 



In order to apply the proposed control design technique, 
let us consider a different basis B such that in the new 
representation p^ = diag ([ l o o o]). This can be achieved by 
considering the Bell-basis 



B 



|00) + |11) |00> — |11) |01) + |10) |01)-|10) 



V2 



V2 



s/2 



V2 



Let B be the unitary matrix realizing the change of basis, i.e. 
pB = pdB. Consider the space decomposition 

77/ = 77s © Jin 

where 77s = span! ^g(|00) + |11))} and Ji R = Jij. We 
have then successfully casted the problem of stabilizing the 
maximally entangled state ( p3| into the problem of achieving 
asymptotic stability of the subspace 77s- Suppose that the 
following generalized measurement is available 

3 

T\p] = Y,M kP Mt 

k=l 

with operators (represented in the computational basis): 



Mi = — (o+ ® J) , M 2 
V4 



1 



(IO(T + ), 



M 3 = \ll - M\M x - M\m 2 



where a + = [§ J], These Kraus operators may be used to 
describe a discrete-time spontaneous emission process, where 
the event associated to M\ 2 corresponds to the decay of one 
qubit (with probability i each), and we neglect the event of 
the two qubits decaying in the same time interval. In the Bell 
basis, the operators take the form 



Aff = 



Mi 



Mo = 






1 

4 

1 

4 



1 1 
4 4 
1 1 
4 4 






0.8536 0.1464 
0.1464 0.8536 
0.8660 

0.8660 



Let us then apply the algorithm developed in section |VT] 
The canonical QR decomposition of the matrices M k returns 
the following triangular factors (we do not report here the 
corresponding orthogonal matrices Qk, see |24} for the final 
form of the controls): 



Ri 



R* = 



' V2 _ s/2 

4 4 





. 

0.8660 0.2887 

0.8165 

0.8660 

0.8660 









V2 




4 


4 















!?2 



4 



V2 

4 





\/2 \/2 
4 4 





According to the proposed approach, by inspection of the 
upper triangular factors i?j we can decide about the feasibility 
of the stabilization task. Indeed, as the blocks Rp t k, k = 
1, . . . , 3 are non-zero blocks, namely 



Rpa 



V5 



Rp,: 



71 



Rp.: 



0.2887 



then the stabilization problem is feasible. 
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Moreover, notice that at this step no further transformation 
is needed on the matrices, as the obtained R factors are already 
decomposed according to 



where — p| fe ker Rp t k- Continuing with the iteration, we 



have then to determine the subspace TL^ = f] k ker Rp' k - By 
inspection one can see that this space is empty, and therefore 
the iteration stops successfully. The set of unitary controls 
that have to be applied when the corresponding outcome k is 
measured is then 



that is: 



' 





V2- 


2 


2 


V2 
2 





V2 
2 





1 








-1 


. 



U k = BQ\B\ 



U 2 



' V2 








v/2" 


2 


2 





1 








V2 
2 








V2 
2 


. 





-1 


. 



0.9856 0.1691 
10 
10 
-0.1691 0.9856 



(24) 



It can be shown by direct computation that the Hamiltonians 
needed to implement these unitary transformation (using ide- 
ally unbounded control pulses in order to make the dissipation 
effect negligible on when the control is acting) form a 3- 
dimensional control algebra |[8l . 

VIII. Conclusions 

Completely positive, trace-preserving maps represent gen- 
eral quantum dynamics for open systems, and if the environ- 
ment is memoryless, also represent generators of discrete-time 
quantum Markov semigroups. Theorem [2] provides a charac- 
terization of the semigroup dynamics that render a certain 
pure state, or the set of states with support on a subspace, 
attractive, by employing LaSalle's Invariance Principle. In 
order to exploit this result for constructive design of stabilizing 
unitary feedback control strategies, we developed a suitable 
linear algebraic tool, which holds some interest per se. We 
proved that a canonical QR decomposition can be derived 
by specializing the well-known orthonormalization approach, 
and that it is key to study the potential of the feedback 



control scheme presented in Section VI In fact, we determined 
which quantum generalized measurements can be simulated 
controlling a given one, and which pure states or subspaces 
can be rendered globally asymptotically stable. Theorem [5] 
gives a constructive procedure to build the controls, and also a 
simple test on the existence of such controls: If the algorithm 
does not stop on the first step, then the control problem has 
a solution. We believe that the provided results also represent 
a mathematical standpoint from which interesting, and more 
challenging, control problems can be tackled, in particular 
when the control choice is constrained by a multipartite 
structure of the system of interest. 

Appendix 

A. Quantum Measurements 

1) Projective Measurements: In quantum mechanics, ob- 
servable quantities are associated to Hermitian operators, 



with their spectrum associated to the possible outcomes. 
Suppose that we are interested in measuring the observable 
C = J^i c i^i- The basic postulates that describe the quantum 
(strong, projective, or von Neumann's) measurements are the 
following: 

(i') The probability of obtaining Cj as the outcome of a 
measure on a system described by the density operator p 
is Pi = Tr^n,). 

(ii') [Liiders's Postulate] Immediately after a measurement 
that gives Cj as an outcome the system state becomes: 

P\* = tracefni/an,) ^iP^i- 

Notice that the spectrum of the observable does not play 
any role in the computation of the probabilities. 

2) Generalized Measurements: If we get information about 
a quantum system by measuring another system which is 
correlated to the former, the projective measurement formalism 
is not enough, but it can be used to derive a more general one. 
A typical procedure to obtain generalized measurements on a 
quantum system of interest is the following: 

• The system of interest A is augmented by adding another 
subsystem B, initially decoupled from A. Let pa® PB, 
with ps = \4>)(4>\, be the joint state; 

• The two systems are coupled through a joint unitary 
evolution Uab\ 

• A direct, von Neumann measurement of an observable 
Xb — ^jUj ) IT, = is performed on B. 

• The conditioned state of the joint system after the mea- 
surement is then of the form 

l -{iA®n j )u AB {p A ® PB )u\ B {i A ®n J ) = P ' Aj 



PAB\ 



Pj 



with pj the probability of obtaining the j-th outcome. 
• One can compute the effect of the measurement on A 
alone, which is nontrivial if Uab entangled the two 
subsystems, i.e. Uab{pa® Pb)U\b cannot be written in 
factorized form. One then gets that p' A ■ = — M^p^Mj, 
with Mj = {£j\UAB\tj>). 
If now the average over the possible outcomes is taken, we 
obtain a state transformation in Rraus form. This construction 
is actually general, in the sense that if the dimension of B 
corresponds (at least) to the necessary number of outcomes, 
any Rraus map can be actually generated this way. 
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